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ABSTRACT 

The present document is a final technical report under the NCC-1-233 research program 
(dated September 15, 1998; see Appendix 5) carried out within co-operation between United 
States' NASA Langley RC and Russia's Goskomoboronprom in aeronautics, and continues 
similar programs, NCCW-73, NCC-1-233 and NCCW 1-233 accomplished in 1996, 1997, and 
1998, respectively. 

The report provides results of "The study of stability of compression-loaded multispan 
composite panels upon failure of elements binding it to panel supports"; these comply with 
requirements established at TsAGI on 24 March 1998 and at NASA on 15 September 1998. 
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INTRODUCTION 

The previous studies in [1, 2] paid major attention, firstly, to postbuckling behavior of 
composite panels and, secondly, to the problem of evaluating strength and buckling of a built-up 
structure with consideration of postbuckling deformation of some components. These studies 
relied on the assumption that the thin-walled elements are attached to supports by links which do 
not fail until the structure becomes destroyed. 

However, the links may turn out to have inadequate strength, and the structure can fail 
due to fracture of the links between the panel and supports. This outcome is evidenced by static 
strength testing of a number of real structures. Figure 0.1 demonstrates airplane wing fracture; 
one can see that the upper, compression-loaded panel has failed because of fracture of links 
between the panel and ribs when the load was much less than the critical value obtained under 
the assumption that the buckling nodal lines rest on the ribs. 

The present research suggests the approach which enables predicting the panel failure 
load as a function of strength of particular links between the panel and supports. 

Let us detail now the essence of the problem. 

One knows that the major load for upper panels in a large-aspect-ratio wing is the 
longitudinal compression caused by general bending. Therefore, parameters of upper stiffened 

panels (including the rib spacing) are specified on the basis of critical stress resultants N* r 
which are computed under the assumption that ribs are nodal lines for the buckled panel. 

This approach is valid if strength of links between the continuous panel and rib/spar 
flanges is sufficient for the buckling mode above. 

However, if not, the panel may buckle at a load N xo less than , upon failure of links 
along one or more supports. 

Thus, there appears the problem of predicting the necessary strength of links between a 
compression-loaded panel and the supports (i.e., ribs and spars) - such that the panel buckling 
nodal lines be on supports; this problem is rather challenging in case a high stress must be 
applied to a monolithic composite stiffened panel. 

The problem is treated here by using the energy method proposed in [3] for analyzing 
behavior of a compression-loaded composite plate with delamination. 

Note that analyzing the strength of links between the panel and intermediate supports is 
different from the usual problem of searching for the necessary stiffness of compliant support; 
the latter formulation proceeds from buckling shapes caused by mutual deformation of the 
supports and permanently attached panel. 


Part 1. The method for studying stability of compression-loaded 
stiffened multispan composite panel after failure of elements 
binding it to transverse and longitudinal supports 

1.1. Solution to the problem of postbuckling behavior of 
unsupported part in compression-loaded multispan panel upon 
failure of support links 

Let us consider a rectangular flat orthotropic panel incorporated in a thin-walled wing 
torsion box (Fig. 1.1). The panel is attached to rib flanges, and its sides, to spar flanges. Panel 
parameters are 

- L , the length, 

- b , the width, 

- a , the rib spacing; 

clearly, L = k m a where k m is the total number of panel bays. At the ends the panel is 
compressed with stress resultants N x0 . 

The subcritical (flat) state of the panel (the state I in Fig. 1.2) is characterized with the 
resultant N xo and the in-plane displacement A x . 

At a certain value of N xo the panel can buckle over a certain length £ = k a (where k is an 
integer), which process makes the panel be separated from k - 1 ribs - refer to state II in Fig. 1 .2. 

The other buckling mode is panel buckling over a certain length £ with separation from 
k - 1 ribs and two longitudinal supports (that is spars) - refer to state IE in Fig. 1 .2. 

If the panel is connected to ribs and spars by rather strong links, then the major buckling 
shape is the usual sinusoidal surface whose nodal lines are on ribs and spars — see state IH in Fig. 1 .2; 

in this case the critical load N CT is calculated by usual formulas. 

Let us address now panel buckling with transition from state I to state II (Fig. 1 .2). 
Assume that the panel in its subcritical condition is compressed with the stress resultant N xo ; 
the panel ends have a mutual longitudinal displacement A x . The panel energy U is equivalent 

1 

to the triangle area — N X0 A X , Fig. 1.3. Upon failure of rib links the panel buckles over a 

portion of length £ , and the load gets decreased to N xo . With the state II implemented, the 
displacement A x does not change. The total initial compression energy U transforms into 
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- compression energy U p and 

- bending energy Ub , 

and the panel has an out-of-plane displacement f . The difference of energy between states I and H, 

R = u-(u p + u b ) O* 1 ) 

is released when the links fail. Energy balance equation (1.1) is represented graphically in Fig. 
1 .3; here, R is the work for destroying the links between the panel and ribs during buckling over 

the length t- ka . 

Thus, energy balance equation (1.1) and the condition A x = const for the transition 

from state I to state II make it possible to determine N xo and the final out-of-plane 
displacement f. The buckling surface length f = ka is evaluated by minimizing N xo with 
respect to k. 

We should relate the compression energy U and the mutual displacement A x to other 


variables in the problem. 

Strains and stress resultants in the midsurface of the orthotropic panel are interrelated as 
follows: 

e x = An N x + Aj2 Ny, £y = A]2 N x + A 22 Fly, Y = A33N X y, (1-2) 

where An, A 22 , A 33 and A 12 are cocos of compliance of the panel in its plane. The mutual 

end displacement may be expressed as, 

A x = L A], N x0 . 0-3) 


Compression potential energy is, 

U = Lb N xo . (1-4) 

To compute the bending energy component Ub and the compression energy component 


Up , we consider the panel buckled over the length i . Assume that the plain (not buckled) parts 

of the panel clamp the transverse edges { x = 0 } and { x = £ } of the separated part. 

In case the panel is tom off the ribs only, the out-of-plane displacement is expressed as, 


2 7tX . 

w = fsin —sin , 
( b 


7iy 


(1.5) 


The midsurface stress resultants are determined by solving the deformation compatibility 


equation from [4]: 
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j 2_ a2 

dy 2 ' dx 2 2 

where we use the nonlinear operator 

' 5 z w 


( 1 . 6 ) 


^(w,w) = 


^ a 2 „, V a 2 w a 2 w 


3x 3yJ 


ax 2 ay 2 


and midsurface strains are 

au i f aw 


\ 2 


8 X = 


• + ' 


dx 2 v dx) 


d\ 1 
e v = — + — 
y ay 2 


aw 

Idy 




au av aw aw 

Y =T" + T _ + T _ ^ — • O- 7 ) 

ay ox ox a y 


Now, stress resultants N x , N v and Nxy are expressed in terms of the function O: 


N x = 


a 2 o 


dy 1 


a 2 ^ a 2 ^ 

Nv — O j Nxy 


,y " ax 2 


dx dy’ 


( 1 . 8 ) 


and substituted (with e x , e y , and y from (1.7) and w from (1.5)) into the formula (1.6) to 
obtain the following equation: 


where 


4(0) = | e /(w,w), 


( \ d 4 d ‘ 4 

A\ ) = A 22TT + 2A 3 


ax" 


ax 2 ay 2 


+ A 


11-4 

dy 


(1.9) 


( 1 . 10 ) 


A 3 =Ai 2 + 2 A 33- 

The relation (1.5) is to be substituted in the right-hand side of (1.6) to obtain 

1 / x f 2 7t 4 2 ( 

-4w,w) = T -j-a ( 


27rx 4 7ix 2 Tty 2 ttx 2rcy 

cos — — - cos — — + cos—: — - cos — — cos 


where 


l 


t 

a = ¥' 


e 


(l.ii) 


The general solution to (1.9) looks like this: 

= ..2 


0 = 0] -N xo 


2 ’ 


( 1 . 12 ) 


where 


r2 T 


<D, = — 

1 32 


2 tux 1 

a 5? I cos— — cos 


4tvC] 8i 

! T-J + ^ 


6 , = 


5, = 


6, = 


2ny 2 2 roc 2 rry 

cos— — -a o 3 cos cos— — 

U *7 l) U 


1 


L 11 


l 22 


3 — ^ A 

A 22 + 2A 3 a +A]|Oc 
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Taking into account (1.12), the stress resultants N x , N y , and in (1.8) can be 
represented in terms of coordinates: 

N x = -N xo - 1 2 (" 6j cos-^ - 5 3 a 4 cos^^j , 


where 


Ny = -t 2 a 2 


t b 

- , 2nX 1 4nX > \ 27CX 2ny 

6 2 1 cos— - - cos — J - 6 3 cos— cos 


XT *2 3s . 2nx . 2ny 
N xy = t a 5 3 sin— ^ sin—, 


(1.13) 


2 f 2 K 2 

1 8 £ 2 


The equation of bending for the orthotropic panel in question may be written on the basis 


of [4]: 


S 4 w 


3 4 w 


5 4 w 


D 11 „ .4 +2D 3 ^2 a. .2 +D 22 a .. 4 a __ 2 a .2 


„ T 3 2 W ^w a 2 w 

-Nv — ~-N„ — r-2Nw - — — = 0, (1.14) 


5x“ J dx 2 5y 2 5y 4 *' x ax 2 ay 2 * dxdy 

where Du , D 2 2 , and D 3 are the panel bending stiffness parameters; these could by using [2] 
be related to stiffnesses of the orthotropic skin and stiffeners. 

The dependence of panel out-of-plane displacement on the external load N xo is 
evaluated by solving Bubnov-Galerkin equation: 

I = J { X sin 2 —■ sin-^- dx = 0 , 

00 to 

where X is the left-hand side expression in (1.14). By allowing for (1.5) and (1.13), we derive, 
-(8Dn+2a D 3 j cos—— sin— + a 4 D 22 sin 2 — sin— 


2 l b 
1 0 0 


n 


t 1 L 


N xo +2t 2 [ 5j -5 3 a 4 cos 2 ^lcos 2 -^ 

£ J b 


2 tux . Try 2 4 
cos — sm - - t a 


( 2 nx 1 4nx 

®2 ^COS ~ — "7 COS 


t 


t 


2kx 2ny 
- o 3 cos — —cos— — 
5 t b 


c ; n 2^Ei e - 2 4s • 2 2tux . 2ny nyl . 2 tux . ny 

sin ^ sin ^ -2t a 8 3 sin ^ sin— cos— |sm 2 — sm — dxdy=0. 


By estimating the integrals and stipulating f * 0 , the following relation appears: 


2 f I ^ x _ 2 

D i i + - a 2 D 3 + — a 4 D 22 J - N xo + y 


n 

£ 2 V 


8l+ a4 (^6 2+ i5 3 ) 


= 0. 


y 


Hereinafter, the new symbols are utilized: 


N n _ ?? p ii , k n k n = 4- 

1N cr - 2 K cr> I '-cr 2 

a k 


1 2 D 3 3 4 D 22 ^ 

1 + — a 3 -^- + — a — 8 — 


V 


D„ 16 d„; 


In these terms, 


a I Di 1 

0 ia - k bM 4 v 

54 = K 5 l+ Ti a452 + I a453 )’ 

t2 = 5:( Nxo " N °)- 


(1.15) 


(1.16) 


The mutual in-plane displacement of the edges { x = 0 } and { x = L } is, 

L f du ’ 

A x =- J — dx. 


dx 


Herein, the component from (1.7) 


du 


____ J_ 

dx ~ H ~ 2 


( dvf^ 2 


V dx) 

should be substituted; and formulas (1.2), (1.5), and (1.13) have to be taken into account to 
produce 

= L = , 9 f 2izy . 9 2ux . 9 ny 

A x = — N x0 + J t 2 I cos ~~~ + 4 sin 2 


t Sin ~b~’ dX> 


or, substituting t from (1.16): 


A x = ^f(5 4 +i)N xo -i N ” 

Oa 


(1.17) 


l_ r (- ' 

>4 

where 84 = 5 4 / 5 1 and £ = t / L . 

To determine external load N xo applied to the buckled portion, we make use of the 
condition A x = const in case the panel comes from state I into state II. Allowing for (1.3) and 

(1.17), the result is, 


N xo - 


5 4 N xo +*N 


n 

cr 


8 4 + ^ 


(1.18) 


Potential energy U of the compression-bent panel is to be subdivided into the bending 
energy component Ub and the midsurface deformation energy component U p : 


U = U b +U p 


(1-19) 
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The bending energy component may be written as 


n j/ f b f L d 2 w] 

Ub ~ 2 J J ° n TT 
z o o [ vax ; 


d 2 wj <3 2 w <9 2 w ^ f* 2 wl _ _ f 5 2 w 

11-2 + 2 ^12 2 -v 2 + ^22 T + ^^33 H — T~ dxdy . 

\dx 2 J dx 2 dy \dy 2 J ydxdy) 


By substituting w from (1.5) into the above expression, we obtain 


u b- 2 4 f 2 J J (D 11 cos 2 ^-D 12 a 2 cos^sin 2 ^ + |D 2 2a 4 sin 4 ^] sin 2 ^ + 

^ o o l v 1 * l 4 l J b 


2 . 2 2 TUX J 7Uy 

+ D 33 a sin — — cos — dx dy . 
£ b 


Upon estimating the integrals and taking into consideration (1.15) and (1.16) the energy is 


^ b= I7 N “( N *»- N ")- 


Or, with (1 . 1 8) in mind, we can write 


^ b _ tt N vn -NE 


Ub = T X N - 


II ^xo ^cr 
CT F 7 

64 + £ 


(1.20) 


The midsurface deformation energy component is represented by 
1 L b 

Up = 2 f | V^l l N x 2 + A 22 N y 2 + 2 Aj 2 N x Ny + A33 N^ 2 j dxdy , 
here, equation (1.13) for stress resultants is substituted: 

U - — - M 2 — ( b f f 5 2 2 *y . f 2 27tx 1 7 4«'| 

Up - 28| N xo + 2 J j 1 81 cos — + 8 2 [cos — + - cos 2-_j + 

+ 5 3 a 4 [(a 4 A, 1 + 2a 2 A 12 + A 22 )cos 2 ^cos 2 ^ + a :! A33Sm 2 ^sin 2 ^]ldxdy. 

l i b t b JJ 

Upon estimating the integrals and taking into consideration (1.16) the energy expression 
becomes 


U n = 


Lb - =2 


P 25 ^ 4 ^ xo + ^ ( U X o “Njj J , 


or, due to (1 . 1 8), the relation ( 1 . 21 ) appears: 


Lb ^4 N xo + £ (nJJ. j 

p " 257 1771 


( 1 . 21 ) 


Now substitute (1.20) and (1.21) in (1.19) to derive the total potential energy for the 
panel in state II: 




( 1 . 22 ) 





11 


Relations (1.4) and (1.22) should be inserted in (1.1) to obtain 

It -J n xo -N°) 2 =R. (1.23) 

We may express work R released during fracture of links between the panel and k - 1 


ribs: 


R = y f b(k-l), (1.24) 

where y f is the specific work necessary to separate the panel from a rib over a unit length of the 
rib flange. 

For the moment when the panel becomes separated from ribs, the stress resultant N x0 is 
determined from (1.23) and (1.24): 


N X o=n5+, 


Ik-1 

2 (64 +f)y f 

1 k 

aA H 


(1.25) 


Both left-hand and right-hand sides in (1.25) could be undimensionalized by dividing 
these by the critical load =N^ for buckling between ribs (refer to state IV in Fig. 1.2); the 
critical load is obtained in the usual way by approximating the out-of-plane displacement with 


, , 7TX . Tty 
the formula w = f sin — sin — : 

a b 


N 1V - ” E il k iv 

^cr — 2 cr > 

a z 


k IV =l + 2a 2 -^- + a 4 ^ 
cr 1 + Za °Dn °D 


11 


11 


where 


a o =“ 


a 

b 


0 < a 0 < 4 4 


D 


ll 


D 2 2 


Upon these transformations, obtain 


N = N n =N° +^p n X, 


where 
- N 


N = 


xo 


N 


IV 

cr 


— n nE 

N — 

iN cr — „ ,jy 

cr 


N 


■ x= iT’ ■^p=I a u( n 5') 2 . p” = 


k - 1 


84 + k ) 

K m ' 


(1.26) 


(1.27) 


(1.28) 


.(1.29) 


The y p variable is the specific energy that would be accumulated in a panel with a unit 


length, which is compressed with the stress resultant N„ ; also, the undimensionalized value X 
is the relative energy spent to separate the panel from a rib. 

In order to determine the minimum possible value N=N n we must minimize (1.28) 
with respect to k , where k = 2, 3, ..., k m . 
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Now we should analyze panel buckling with transition from state I to state HI (Fig. 1 .2) - 
that is, panel buckling with separation over the length £ = k a from k - 1 ribs and two spars to 
which the panel is attached by its side edges. 

Upon the separation from spars the longitudinal edges are free. Displacement W may be 
expressed as, 


f . 2 ** 
w = f sm — . 


(1.30) 


With links between the panel and spars broken, the panel buckles over the length £ , and 
load decreases to N xo . During transition from state I to state III the displacement A x does not 
vary. So the panel buckling critical stress resultant is 


n iii = ZL.. d ii , m . m _ J_ 

1 cr 2 ^cr » 1 '-cr — , 2 " 

a k 


(1.31) 


Now we use analogy with theories of [5] for the angle of the tangent to the midsurface of 
the buckled panel in the longitudinal direction and take into account (1.30) to obtain 

S = 0 8, 

where 


cr 9 


. Nxo-nJ? 

0 = J 8 CT 


N 


m 

cr 


0 CT = sin 


271X 


To derive the displacement w , we utilize the boundary condition {w = 0 at x = 0 and x =£ }: 

= 0 fs CT ds = 0 — sin 2 — . 

J CT TT o 


w 


This relation is compared with (1 .30) to have, 


^fli 2 _Nxo-N® 


8 £ 2 


N 


m 

cr 


(1.32) 


The stress resultants N x , N y , and N^ for state III are 


N x =-N 


xo> 


N v = N vv = 0. 


'y - ‘’•xy ~ (1-33) 

The mutual displacement of edges { x = 0 } and { x = L } is written with due 
consideration of (1.2), (1.5), (1.33), and (1.32): 


t h , n xo -n!? 

A x = -J— dx = 2^ xo CT 

ex 


N 


m 

cr 


+ LAuN X0 . 


(1.34) 


13 


The external load N xo of the buckled panel may be evaluated by using the condition 
A x = const for the panel transition from state I to state in. With (1.3) and (1.34) in force, we 
have 


Nxo — 


where 


From (1.32) and (1.35): 


t 2 = 


^ xo T ^cr 

e a +2t 

(1-35) 

*cr • 

(1.36) 

Nxo - N™ 
*1 ( e cr +2l) 

. (1.37) 


The bending energy component is expressed as. 


= b \ 

Ub = — jD n 
z 0 


W 2 


vax 2 ; 


dx. 


Now substitute w from (1.30) and allow for (1.36) and (1.37) to obtain 

tin 

^XO ~ A1 

*cr 


u b = ^bAr L 4NF Nxo N * 


e cr +2^ 


(1.38) 


The relation (1.33) is now taken into account to derive the compression energy 
component for the plate midsurface: 

b h y 
Up = T K,N 2 dx. 


Herein, we substitute stress resultants from (1.33) and (1.35): 


U p = Lb^- 


E cr+ 2 f 


(1.39) 


With relations (1.4), (1.38), and (1.39) obtained, equation (1.1) may be re-written as 


£ b 


L ll 


e cr + ^ 


(n X 0 “ N™ ) = R 


2 ( e <*/2 + ef 

Real panels are characterized with the interrelation e cr «£ , so 


-cr 

\2 


Lb^-(N X0 -NS') 2 = R. 


(1.40) 


Let us decompose R into two summands: 


R — Rf + R r , 


(1.41) 
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where 


Rf = Y f b(k-l), R r =2 y r ak, (1.42) 

y r is the specific energy necessary to separate the panel from a unit length of the spar. 

Here, we substitute (1.41) and ( 1 .42) in ( 1 .40) and take into analysis that i = k / k m : 


N “ = N “ + ] (nr Tf +2o » r ') ■ <‘- 43) 

Both left-hand and right-hand sides of (1.43) may be divided by the critical stress 
resultant N CT =N™ for panel buckling within a single bay with no separation from ribs and 
spars (refer to state IV in Fig. 1.2): 

N = N m = N? + p m X , (1.44) 

where X corresponds to (1.29) and 


X r = 


y? 


N ffl 

iN cr 


N 

N 


m 

cr 

IV 

cr 


P m = 


~[k(l + 2<x„X r )-l], X r = 

L m 


^ 2k 
x 


(1.45) 


For particular values of a 0 and ratio X r we should determine the value of k (with 

k = 1, 2, ..., k m ) which provides the minimum N level. 

The maximum load N that may be carried by the panel not separating from supports 
(with panel parameters and link fracture energy prescribed) corresponds to the lower value 

among N n and N m (as provided by (1.28) and (1.44)) while assuming that N < 1 . Note that 
both (1.28) and (1.44) have been obtained under the assumption of N xo < ; so the relations 

are allowed to be united as 


N 11 

if 

n" £N m and N x0 < n" , 

N m 

if 

N'^N 111 and N^sN^, 

1 

if 

. (r;n rrlll) ", 
nun IN ,N 1 >1 . 


(1.46) 


Let us use X CT to denote the minimum possible value of the relative link fracture energy 
X for links between the panel and a rib; this X CT value at the X r / X ratio specified ensures 

N = 1 (see (1.46)), and no greater X value could increase the failure load. From (1.28), (1.44), 
and (1.46) we obtain 


'x? r , if x” > X™ 
'.x™, if xS<x»'’ 
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where 


yll — 
A cr — 


max 

k=2, 3, .... k, 


I-(i-nS) 2 }, {^(i-nS 1 ) 


These formulas utilize the k values which ensure < 1 and ^ 1 , respectively. 

When X > X CT , the panel is attached to supports strongly enough for the panel to buckle 
between ribs (i.e., to come to the state IV) with no separation from ribs and/or spars. Thus, the 
critical value X CT of link fracture relative energy at X r prescribed is an important indicator of 
strength of joints between the panel and its longitudinal and transverse supports. 

In particular, the present methodology may be employed to evaluate link fracture energy 
for a skin and the longitudinal and transverse stiffeners. 


1.2. Dependence of ultimate load on panel parameters and link 
fracture energy available: A parametric study 

In the general case an orthotropic panel may be described as a composite skin (with a 
thickness h and a symmetric stacking) stiffened with stringers. The material characteristics of 

the orthotropic skin include Ef , E§, G°, pf , and - that is, elastic and shear moduli and 
Poisson's ratio. Stringer parameters include the following set: 

- b s = b / (n s + 1) , the stiffener spacing, 

- n s , the total number of stiffeners, 

- E s , the stiffener Young modulus, 

- F s , the stiffener cross-sectional area, 

- I s , the stiffener cross section inertia moment, and „ 

- z s , the distance from the stiffener cross section centroid to the skin midsurface. 


For the panel the compliance cocos are 


_ 1 . 1 
1 1 ~ Eih ’ 


An = — . A 22 =— a.___JLL_a _ ^2 1 


a 12 = -^ = a 21 =-^, a 33 = ^. 


where 


1 + r 


E, =Ef(l + r). E 2 =E|^i, = 

G = G°, ;.r(l-nfnj), r = 


E S F S 

Efhb s 
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The panel bending stiffness is determined through skin and stringer stiffnesses as 
follows: 


D il = D ?i +D n> 


D 22 = D° 


where 


D?, = 


Ef h 3 


11 12(1 


D°22 = 


E 2 h 3 


22 » 


D3 - D3 - D 12 +2D33, 


D33 = Dr* = 


33 


G° h 3 
12 


12(1 - n? mS)’ 

° 11= 1 0 + b S+ b S ( Zs_hl ) ’ D 12 =D ?2 = + M? D 22)> 

1 — Ui Ut u S u S 9 


E?h l2 E«I« E B F e 

hf + 

Mi M2 “s u s 

and hj is the distance from the skin midsurface to the centroid of the cross section of the 
stiffener with the neighbouring skin. 


By analyzing formulas ( 1 . 28 ) and ( 1 . 44 ) it can be concluded that, with X and X r /X 


specified, the value of N (and the integer k which minimizes N) depends on the following 
parameters: 


A 22 A 3 Du Dj 

A 11 ’ A ll ’ D 22 ’ D 


It is easy to write similar parameters for a plain (unstiffened) orthotropic panel: 

a k M o 

b l b l 


In case the panel is isotropic the N value does only depend on a 0 , k m , and p . 

Thus, the total number of parameters influencing N is rather large. In this connection we 
limit ourselves to analysis of a few typical versions that help nevertheless to detail the entire 
methodology. 

Firstly, let us consider a 10 -bay panel (with k m = 10 ) attached to ribs and spars; assume 
that the rib spacing is much less than the distance between spars (here, a 0 = 0 . 2 ). At the 
initiatory stage, we presume in addition that the relative strength of panel/spar links is notably 
greater than that of panel/rib links (X r X » 1 ) — this means that, as the load is growing, the 
panel may separate from ribs only. We will address three versions of panel design: 


- isotropic panel, 

- homogeneous orthotropic panel, and 

- inhomogeneous orthotropic panel (stiffened in the longitudinal direction). 
Figure 1 .4 represents the dependence of N on X for the isotropic panel in which 


17 


D 11 A 22 _ D 11 A 3 

D 2 2 A n D 3 A u 

If the available link fracture energy X grows, the ultimate load N gradually increases, which is 
accompanied with shortening of the separation zone length t = k a . For example, within the 
interval {X 56 <X<X 45 } the panel buckles and separates from four ribs (k = 5), whereas 
within the interval {X 56 <X<X CT }, from three ribs (k = 4). When X^X CT the panel 
buckles between ribs and obtains the shape corresponding to state IV in Fig. 1 .2. From Fig. 1 .4 it 
is seen that the isotropic panel is characterized with X CT = 0.53 . 

Figure 1.5 depicts the same kind of dependence for the homogeneous orthotropic panel 
(that has a composite skin with a symmetric layup and no stiffeners) in which 


D ll _ A 22 
D22 Ajj 


=2.87, 


D 


D 


— =3.64, 



The panel dimensions are the same as in the previous example. Here, we have the critical 
value X CT =0.71, which is higher than that in the isotropic structure. 

Figure 1.6 shows the «N-X» diagram for the orthotropic panel (with a composite 
skin and stiffeners) in which the in-plane dimensions are identical to those in previous 
versions; here, 


D 


11 _ 


D-- 


70.1, 


L 22 


=3.46, 


D 


ll _ 


'22 A 11 D 3 A n 

Here, the structure has the critical value X CT = 0.77, which is slightly higher than that in 
the unstiffened composite panel. 

By comparing values of X CT for the three versions, it is clear that increasing the 
longitudinal stiffness results in greater X CT values. 

Of interest is the dependence of X CT on the relative rib spacing a 0 . 

Figure 1.7 demonstrates how X CT of the isotropic panel depends on a 0 ; Figures 1.8 

and 1 .9 show similar functions for the homogeneous composite skin and the orthotropic panel 
with the stiffnesses reported above. 

From Figures 1.7 through 1.9 it is seen that X CT notably decreases as a 0 grows. 

Below, we consider the same panels but assume that these could buckle in accordance 
with models II and III. The ratio X r /X = Y r /yf for these examples reaches various levels: 

- 6.0 for the isotropic unstiffened panel, 

- 4.0 for the composite unstiffened panel, and 


= 88 . 8 , 


=5.63. 
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- 3.0 for the composite stiffened panel. 

Figures 1.10 through 1.12 represent relevant «N-X» diagrams. The curve corresponding 
to state H intersects (at a certain point X = X m _» n ) with the curve for state m. For the panel 

versions under consideration the X CT value are on the curve corresponding to state II. 

Figures 1.13 through 1.15 demonstrate the X CT dependence on a 0 . As for the isotropic 
panel the curve corresponding to state II intersects the state IE curve at the point with a c =0.23. 
Consequently, if the available X value is below X CT , whereas a 0 <0.23, then the buckled 
panel may separate from ribs only - state II. If, however, the available X value is below X CT , 
but a 0 > 0.23, the buckled panel separates from ribs and spars - state m. And the structure with 
X > X CT buckles between ribs with no failure of links to both the ribs and spars - state IV. As 
for homogeneous composite and structurally orthotropic panels (with the respective X r / X 
values written above) the coordmate a 0 of intersection of curves corresponding to states II and 

HI is slightly larger than that for the isotropic panel (here, a o =0.31 and a o =0.34, 
respectively). 

To conclude the investigation, we carried out studies on dependence of N on X and 
X ct on a 0 for panels with different lengths — at k m = 7, 8, 9, 10. 

Figure 1.16 depicts the dependence of N on X for the isotropic panel at a 0 =0.2. The 
bright circles on the curves denote panel transition from type m buckling to type H buckling. 

Figure 1.17 demonstrates how X CT for the same panel depends on ct 0 . These curves 
suggest that decreasmg the total number of bays in the generic panel decreases X CT , other things 
being equal. This is associated with the fact that, with the total number of bays decreased, the 
potential energy of the unbuckled panel portion (whose length is L-f) decreases as compared 
with the total energy of the entire panel upon buckling. 

Figure 1.18 shows dependence of X CT on a 0 for the composite stiffened panel whose 
stiffnesses are represented in the Figure. 

Figure 1.19 provides the computed results on X CT as a function of the "relative stiffiiess 
capability" in the longitudinal direction - D„ D 22 ; the ct 0 parameter assumes certain values: 
0.2, 0.25, 0.3, 0.35, and 0.4. In the case of D n D 22 >20 the parameter is seen not to influence 
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Part 2. Fracture energy and strength of links between the multispan 
panel and supports: Critical values for ensuring traditional buckling with 
nodal lines being on the supports 

Part 1 above introduced the notion of "specific link fracture energy X " for a compression- 
loaded multispan panel attached to supports; also, Part 1 provides analytical expressions and some 
parametric analyses that establish dependence of the panel failure load on X . To carry out real 
analyses, one should be able to determine X for particular supports and joints. There exist 
various designer solutions in respect of supporting elements and types of panel-to-support joints. 
In addition, the "panel-joint-support" system fails at its weakest point. So the link fracture 
energy is recommended to be determined by resorting to special experiments with real structure 
parts; experimenters should record 

- the load-displacement diagram, «q - A », and 

- the ultimate load q p 


(here, q is the linear density of the peeling-off load normal to the panel surface and applied to 
the line of joint between the panel and the support, and A is the displacement in the q 
direction). This diagram must be obtained while increasing the load q from zero to "panel- 

a p 

support" link fracture. The area outlined by this diagram, JqdA is the equivalent to the link 

0 

fracture energy; here, A p is the displacement at fracture. 


However, the experimenting with real structural parts is difficult to implement for a 
number of reasons. In this connection we propose the method for theoretically determining the 
link fracture energy while assuming that the weakest point in the "panel-support" system is the 
joint (i.e., a bolted or riveted joint or an adhesive bond) whereas the support itself is rather stiff. 
If the support is a beam-like rib transferring its load to spars, then the rib deformation energy, 


U 


b.w. 




aQ 2 (q 

2GF 



( 2 . 1 ) 


can be added to the fracture energy of the joint, as this energy is released at the failure time 


instant. In the above equation the symbols M(q p ) and Q(q p ) are the bending moment and the 


shear force at the load q p , respectively, and EJ and GF are the rib bending and shear 


stiffnesses, respectively. 

Let us preliminarily address a simplified approach to bolt/rivet fracture energy, provided 
that the joint fails because of breakage of the rod of a bolt or rivet. 
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2.1. The method to evaluate link fracture energy for bolted and riveted 
joints. Relation to load-displacement diagram and bolt/rivet fracture 
load 

To carry out real analyses, designers should evaluate parameters y f and y r . We suggest 
the method for determining these on the basis of experimentally obtained «q - A » diagram; 
here, q is the linear density of the tensile load that tears the panel off the rib/spar, and A is a 
characteristic displacement in the joint subjected to tension. The specific fracture energy (per 
unit length of the joint) is, 

a p 

Y=JqdA, (2.2) 

0 

where A p is the characteristic displacement at fracture time instant. 

If the panel is attached to the supports by using bolts and/or rivets, and the fasteners fail 
due to rod breakage, then y may be approximately computed on the basis of the «o-e» 
diagram for the bolt/rivet material. 

Assume that the « a - e » diagram for tension of a bolt/rivet has the shape represented in 
Fig. 2.1; here, we employ the following notation: a 0 2 is the yield stress; 82 = a 02 E ; and a p 

and e p are the ultimate stress and the relative elongation at fracture, respectively. 

With this, the rivet/bolt fracture energy and the y value may be calculated as, 

A = HFS d and y=— HFS d , (2.3) 

c 

where 

H and F are rivet/bolt rod length and cross-sectional area, respectively, 
c is the rivet/bolt spacing in a row, 
n is the number of rivet/bolt rows in the joint, and 

e p 

S d = Jade is the area occupied by the « a - e » diagram; the area may be approximately 
0 

computed as, 

Sd=°pE p Tl, (2.4) 

where r| is the «ct - 8 » diagram filling ratio: 




, 0.5 < r\ < 1 . 


(2.5) 
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Thus, to determine y , we should substitute the integral (2.4) into the formula (2.3): 

y = nq P A p , (2.6) 

note that 

a P F 

q p =-^-n and A p =He p . (2.7) 

By substituting (2.6) into (1.29) and (1.45), we obtain relative link fracture energy values 
X (for ribs) and X r (for spars): 

X _ 2 TlfCi P f A P f X - 2 T 1 r( 1 P r A P f n ox 

A ( iv\ 2 ’ Ar ” / IV \2 ’ ( 2 - 8 ) 

Ah(nS') A n (N^ v ) 

Here, q p f, q p r» "H f , "H r » ^ pf > A F are the characteristics introduced above for evaluating 

joints between the panel and ribs and spars, respectively. 

The values X and X r enable one to use either plots in Figs 1.10 through 1.12 or 
immediate computation in order to establish 

- the load at failure of the compression-loaded multispan panel and 

- the real buckling mode (that is, II, HI, or IV). 

There exists an opposite problem: one may want to ensure such a buckling surface shape 
between ribs that links do not fail; in this situation the condition X > X CT makes it possible to 
determine the required parameters of bolts and rivets, as (2.7) includes fracture characteristics 
a p and £ p , bolt/rivet section area F , spacing c , and number of rows n . 


2.2. Finite element analysis of link fracture energy for bolted joints 

A joint is a system of discrete, irregular components which are in contact. A joint may 

* 

fail not only because of the bolt/rivet rod breakage considered in 2.1 but also due to head 
separation and/or breakage of a plate/support near the bolt/rivet head. In this connection the 
discrete link fracture energy should be established by utilizing Finite Element Methods (FEM). 
The joint with the neighbouring portions of a rib and/or spar should be "dissociated" into finite 
elements, and for each element we must determine the stress tensor cry and the strain tensor Ey ; 

the joint is assumed to be loaded with tension due to detachment of the parts. In this case the 
energy accumulated by the components during the system deformation to failure is described by 
the following expression: 
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C P 

^ ~ j !° ij dsjj dV , .. 

V 0 K ’ 

where V is the total volume of all elements in the unit, and e p is the breaking strain (let us 

assume that the structural system is destroyed if the breaking strain is attained in at least one point). 

It is obvious that the problem of establishing stresses and strains within the unit is 
notably nonlinear - due to two circumstances: 

— the components are in a variable contact and 

- materials are physically nonlinear, which is of importance in structures with heavy 

loads. 

Each of these problems is traditionally solved by employing stepwise approximation 

methods, which require much more computation costs than dealing with linear problems. Here, 

we propose an efficient means for determining contact stresses and allowing for material 
nonlinearity within FEM. 


Method for solving contact problems 

The principal set of equations is obtained by minimizing the Lagrange energy functional 
and looks like this: 


[K]{5} {R}, (2.10) 

Where {8} is the nodal displacement vector, {R} is the external load vector, and [K] is the 

trtructure stiffness matrix. To solve the set (2.10) we must specify boundtuy displacement 
conditions for a certain part of the domain: 


5 = U 0 


( 2 . 11 ) 


Let us consider two bodies of arbitrary shapes (i and'j am indices to identify the relevant 
bodies) in the Cartesian coordinate system OXY (Fig. 2.2). 

Let S lc and S jc be the assumed surfaces of contact, i.e„ those portions of the bodies i 


and) that can interact during loading. The surface S ic at ever, point C* (k is a point number, 
k=l, 2, ...) has an outward normal n* . Points C* and C jk having identical second indices 

will hereafter be referred to as conjugate points if these could establish mutual contact. As noted 
in [6], whether points are conjugate can only be detemtined with a fair degree of accuracy and 
prior to solving the problem if the pattern of contact deformation is obvious and the bodies have 
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rather simple geometries. If it is not the case, the contact area will be outlined by employing an 
algorithm based on criteria of mutual nonpenetration of bodies. 

Introduce radius vectors r oi and r oj to identify initial positions of the points C* and 

C jk , respectively (Fig. 2.2). After loading, positions of these points in the three-dimensional 
space will be defined by the relations, 

{ r ik}={r 0 i} + {6 ik }, ( r jk } ={r 0 j} +{5 jk } , (2.12) 

where {§*} and {8 jk ] are displacement vectors of the conjugate points of the i-th and j-th 
bodies. The criterion of contact of the points C* and C jk can be written as, 


({ r ik} ~ { r jk }) { n ik} = 0 


^ik e ^ic 


(2.13) 

After considenng equation (2.12), from equation (2.13) we obtain the compatibility 
condition: 


(f r &} { r jk}) { n ik} -({S jk } -{5ik}) {njk ]Jc jk e s . c . 

For conjugate noncontacting points the following condition should hold: 


({ r ik} ~ { r jk }) { n ik ) - 0 


C ik eS ic ’ 


(2.14) 


(2.15) 


which, in its essence, expresses the condition of mutual nonpenetration of bodies. In the 
projection on the normal n^ , the condition (2.14) becomes 

8 (n) ,(n) (n) 

*>ik +6 0 -6 jk , (2.J6) 

where 5* and 8 { jk } are displacements of conjugate points of the bodies i and j along the 

outward normal; sj, } is the initial (positive or negative) interference between the conjugate 
points as measured along the normal . The equations of equilibrium of isolated bodies i and j 
(with no contact between these) in block matrix notation are ■ 


"K in 

K i2 i 


K 


i!2 


K; 


i22 


r«ii 

l5 i2 


LQi 


"Kj„ 

L K J2I 


Kj 12 


f5j, 


K 


> = ^ 


J22 




r j 

Qj 


(2.17) 


where 


S n and Sj, are displacements of nodes (on the bodies i and j, respectively), that do not 
contact mutually, 

5 i2 and 8 j 2 are displacements of nodes which contact one the other, 
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Pj and Pj are prescribed external loads, and 


Qi and Qj are contact forces to be found. 


Now we resort to the Gaussian elimination procedure for block matrices to describe 
unknown displacements of contact nodes: 

[K-]{ 6i } = {<}, [K*]{Sj} ={ R *}. (2.18) 

where 

] =[Kj 22 -K^iKj^jKip j, j | = j 2 2 - K j2iK j, l jK j 12 j , 

(2.19) 

[R-j^Qi-KjuKu'.P,], [R j ] = [Q j -K j^K j/jPj]. 

Assume that the conjugate (k-th) nodes of the i-th and j-th bodies are interrelated in 


a local coordinate system: 



where 



[c] -[c ]1 
L-[C] [c]J’ 


( 2 . 20 ) 


( 2 . 21 ) 


Qik and Qj^ are the contact forces applied to the conjugate points, 


[C] = 


[C] 


0 C 


0 




in two-dimensional problems, 


C n 0 0 

0 0 

0 0 c n 


in three-dimensional problems, and 


( 2 . 22 ) 


C n , C ^ , and C ^ are link stiffnesses along the local coordinate system axes n, £, , and 


r| , respectively. 

Let us introduce [X] , the direction cosine matrix relating the local coordinate system 
(n, r|) to the global one (x, y, z). Expressions (2.18) and (2.20) should be summed up to 

produce the equilibrium equations for the bodies in contact: 


i[K?]+[y 0 ] 

1 

o' ’ 

> - J 

1 

r»,i 

r‘ 

1 

1 

o 

‘ ‘ 

[K*] + [y°] 

i 5 )} 

r 1 

L j . 


(2.23) 
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Values in the |y 0 j matrix characterize link stiffnesses in the general coordinate system: 

[Ykij] = [^kij] T [C] [Xfcij]- (2.24) 

If the i-th body has interference, {5 0 } , the right-hand side of the first equation in (2.19) 
should be complemented with the load summand 

{R 0 > =[Ki]{5 0 }. (2.25) 

In case the bodies are not in contact, we obtain |y°J = 0, and the formula (2.23) is re- 
written as, 


K 


M 


Si 

l 5 JJ 


_K i 12 *fn p i 


-K 


-1 


J12 K jn r j 


(2.26) 


The latter system for unknowns {5,} and {5j} may be solved immediately. But when there 
occurs a contact, the system (2.23) is to be solved using the following algorithm 


1 . Specify the contact nodes and their stiffness matrices [c^ ] . 

2. Compute the direction cosine matrix [X] for the contact nodes. 

3. Calculate the substructure stiffness matrices and right-hand side vectors - in 
accordance with (2. 1 8). 

4. Use (2.18) to derive equilibrium equations (2.19) for contact points, by utilizing the 
Gaussian elimination procedure. 

5. Use (2.17) to determine along-the-normal interference 5 R between nodes: 

§Rk=Sjk-S!L-5ok 

6. Find the forces R applied to the body i due to the interference 5 Rk : 
(r j } = [K j ]{5 Rk } ; add the forces to the external load vector of the body. 

7. Utilize expressions (2.24) in order to transform link stiffness matrices into the global 
coordinate system; add these to the structural stiffness matrix (2.23). 

8. From the relation (2.23) determine the contact node displacement vectors (Sj } and 



9. Use (2.22) to compute contact force components along normals: 

Q k =C k X (5 ik -6 jk ). (2.27) 
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10. Check up the nonpenetration condition (2.16) being met. If yes, go to item 1 1 below; 
otherwise, repeat the analysis starting from item 5. 

1 1 . Terminate the iterations and calculate stresses by employing usual FEM relations. 


Allowance for materials plasticity 

If at least one component in the joint behaves nonlinearly, the contact problem is much 
more difficult to solve. To describe physical nonlinearity, we use the strain theory of [7]. 

The problem is solved under the assumption that the structural displacement are small in 
comparison with overall dimensions. To outline the contact area, the above algorithm is 

involved. The coefficients K * of the equilibrium system (2.23) for the interacting elements: 

[Kj(5)]+[y°] — [y °] 

’ -cy°] [Kj(5)] + [y°] 

become dependent on nodal displacements; so these should be determined by resorting to a 
successive approximation procedure in couple with the stepwise (incremental) loading method. 

The algorithm for treating the elastic-plastic contact problem has been detailed in [8]. To 
improve iterations for outlining the plastic deformation zones, we utilize the secant modulus E s 

correction approach similar to the Reference [9] overrelaxation method for linear algebraic 
systems. The correction is implemented as follows: 

^si — ^si "*”®(^si+l — ^si)> 

where i is the plasticity iteration number and co is the relaxation factor which is varied wit hin 
the iterative cycle: 

co =1 + Aco /i a . (2.29) 

To define an optimum co value, additional numerical experiments have been mounted 
for various problems. The parameter aco involved with (2.29) was changed from 0 to 0.8, 
whereas a was specified to be from 0.5 to 4.0. Iterations were terminated when a relative error 
had been less than 0.005. 

It turned out that the minimum number of iterations is required in most problems when 
aco =0.6 and a=2.0. These values were utilized to solve the problems reported here. 



r*.i 
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(2.28) 
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Example problems 

To improve (and approve) the present method for solving contort A f 
we analyzed stresses and strains in the bolted join, model to, was studied iTT' 0 " 
optical polarization devices, as described in [10). A schematic of to boTITto I" d 

application may be seen in Fig. 2.3. * d ™ 0ad 

The flat models of bolts were manufactured out of 4 2 mm 1 , i 
active material. OASO. The bolt model 1 was placed on a ° P ' iCal,y 

bar 3. The right-hand side of L F T ofb °'< '• support 2, and 

force P had 1 appfild $ — « - - head afier to 

« are r ot: ":ri r: z r ,oad ~ * — *> - 

,'ctoes are easy compare to co'nclnde to £ ^ 

converge well in both the pattern and amplitudes of maximum shear stresses ? ^ 

F^ne 2.6 represents theoret , a| variation of ^ s(ress coMon 

o " -P/2 JT A ^ B ' dePeDdi ” 8 °" - ^ *** M, here, 

non, /(td) and a, ,s to principal tensile stress at a point In the F 

experimental data for h d - 0.4, 0.5, and 0.6. The most severe s, ^ 

^ fillet benveen the bol, rod and head Th f • ^ ^ POin, B " “ 

should be expected to occur ; a, a higher load to Z mL""' ^ 

-ed rt: on a “ — 

ouuu kg through two bolts The FF t r . 
depicted in Fig. 2 7b- it inclnd,. . ' E del of ,he structure is 

g 4,/b, ,t includes simplex axisymmetric finite elements The w„ u 

unknown displacements is 1298. When solving die cento nroble 

displacements over to zone where to boh interacts with 

Computed results are compared with experimental data obtained^ st^ah^ga^ne ^ ^ 

Figure 2.7c demonstrates axial stress profiles over to external H ' eS ‘ S ' 

the cylinder. Points show the experimentally obtained values The theory ^ 
regarded as converging well. ™ * d ex P enmen t may be 
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Computing the link fracture energy 

The link fracture energy can be determined by using (2.9). Integrals in (2.9) may be 
estimated by the following procedure. The external load in subdivided into a number of 
increments. As the load grows, each element in the unit is evaluated in what concerns stress and 
strain tensor components; their algebraic products are summed up for each increment by 
utilizing the trapezoidal technique. The volume integral is the sum over all element of the joint. 

Let us address tension of an aluminum bolt whose dimensions are represented in Fig. 2.8. 
The bolt material has the following characteristics: E=7200 kg/sq.mm, ultimate stress 
Op =38 kg/sq.mm, assumed yield limit gq 2=24 kg/sq.mm, and elongation at rupture, e p =0.1. 

The tensile «(7 — s» diagram of the material (see Fig. 2.1) is modeled by two straight 
lines: the first line approximates elastic bolt behavior, whereas the second line runs from the 
assumed yield limit point to the ultimate state point to approximate behavior of materials in the 
case of notable plastic deformation. 

Figure 2.8 depicts dependence of bolt fracture energy on stress a . We represent three 
versions. In the first version the bolt head rests on the support (b) manufactured of the same 
material. In the second version the head is in contact with an absolutely rigid support (c). In the 
third case we consider the bolt rod only. The analytical model for the joint under consideration is 
similar to the one demonstrated in Fig. 2.4 - except for the fact that, instead of flat (two- 
dimensional) finite elements, we used axisymmetric ones. Established as a failure criterion is the 
strain intensity e p =0.1. In Fig. 2.8 it is seen that the link fracture energy (symbolized with 

crossed squares) is slightly greater for the bolt contacting with the deformable support ( 1 ), than 
the value for the bolt resting on the rigid support (2). However, link fracture energy for the two 
versions is much less than that for the separated bolt rod (3); this is so due to the considerable 
stress/strain concentration at the bolt head-to-rod fillet. 

Now address the second example: a skin 2 attached to a rib or spar 3 by a single-row 
screw joint; one of the screws 1 is depicted in Fig. 2.9. The skin thickness is 6 mm, the screw 
diameter is 8 mm, and the screw rod length is 20 mm. 

When computing, we assumed that the rib/spar 3 is made of aluminum alloy, the screw 
is made of either steel or aluminum alloy; the steel has E=21.000 kg/sq.mm, c p = 60 kg/sq.mm, 

CT 02 = 40 kg/sq.mm, and e p =0.1. The skin is made of aluminum alloy in the first version, 

and of a composite material in the second version. The composite skin includes 24 layers each 
0.25 mm thick. A substack of 4 layers has a quasi-isotropic layup (see [11]): [0, +45, -45, 90], 



Characteristics of a layer may be seen in the mid-column of Table 2.1; the subscript 1 
designates the fiber direction, and the subscript 2, the transverse direction; p stands for tensile 
load, and c is for compression. Characteristics of the entire stack are reported in the rightmost 
column of Table 2.1. 


Table 2.1 



Layer data 
kg/sq.mm 

Stack data 
kg/sq.mm 

El 

13000 

5200 

e 2 

700 

5200 

G 12 

590 

2000 

Pl2 

0.36 

0.29 

°lp 

153 

51 

a 2p 

3.4 

51 

CT lc 

120 

40 

ct 2c 

3.4 

40 

X 12S 

6.0 

30 


The metallic materials (that is, steel and aluminum alloy) are assumed to become 

broken when the strain reaches its limit; for the composite material (whose fibers and matrix 

behave linearly to failure) we use the Tsai criterion from [12]: 

2 2 2 
Gi _ q 1 q 2 02 *12 . 

S? si sj S 2 

where S = Tj 2 s is the limiting shear stress, and Sj and S 2 are the ultimate normal stresses along 

* 

and transverse to fibers, respectively: 


Si 


a lp at CTj >0, 
ai c at aj <0, 


S 2 = 


ct 2p at a 2 >0, 
cj 2c at ct 2 <0. 


Mechanical properties of the quasi-isotropic plate are obtained by using the Tsai method in [13]. 

Figure 2.10 demonstrates the deformed model (composed of axisymmetric finite 
elements). Figure 2.11 represents dependence of link fracture energy on the force P (which is 
calculated by integrating the distributed load p over the skin surface area). 




From Fig. 2.11 it is clear that link fracture energy for the aluminum skin and the steel 
screw (b) is approximately twice that for the structure with an aluminum screw (a). Failure 
mode implemented is also depending on the screw material: the aluminum screw fails due to 
head separation, whereas the steel screw survives the skin. 

The composite skin (represented by lines c and d) gets broken at the screw head when the 
separation force is at a rather low level. The fracture energy in this structure may be increased to 
the line b by resorting to usual designing solutions (detailed, for example, in [14]): increasing 
the screw head diameter, increasing the plate thickness around the joint area, etc. 


2.3. Establishing parameters necessary to ensure a specified 
panel-to-support link strength 

Consider panel separation from support in case the support (rib) deformation energy is 
comparable with the link fracture energy for a bolted or riveted joint. The work for separating 
the panel from a rib may be written as, 

Yf=Y P f+Ybf> (2.30) 

where y p f is the specific work of breaking the fasteners between the panel and the rib (the work 

is referred to a unit length of the rib flange) and y bf is the specific potential energy of rib 
deformation at the panel separation instant. 

We can model the rib as a beam simply supported by spars at edges {y=0} and {y=b}. 
At the moment which immediately precedes the panel separation from the rib, the latter is 
loaded with a uniformly distributed transverse force q p . The rib deformation potential energy 

Ubf is computed by using (2.1) where 


M = -^-(y 2 -yb), Q = (2y - b) , , (2.31) 

a f is the coefficient depending on the rib cross-section shape; the coefficient may be assumed 
to be equal to 1 .0 for a web-like rib. 

By substituting (2.31) into (2.1) and integrating, we obtain the following expression that 
relates the specific potential energy Y bf = Ubf /b to the load q p : 


where 


Y bf = c bfq p 


c bf = 


240(EJ) f 


1 + 


10 (EJ) f 

b 2 (OFw) f J 


(2.32) 

(2.33) 



The coefficient c^f quantifies the rib compliance under in-plane bending and shear. 

The specific fracture work (referred to a unit length of the joint) can be written as 
follows: 

r P f=c pf q*. < 2 - 34 > 

The coefficient c pf quantifies the compliance of fasteners between the panel and the rib when 
loaded to failure. A particular expression for c p f depends on both the panel/rib joint type and 

the failure mode. For example, if a panel is attached to supports by bolts or rivets and if a bolt or 
rivet is destroyed due to breakage of the rod, then, in accordance with (2.6) and (2.7), we obtain 

nfCfHf (2.35) 


Cpf HfFpfEspf’ 

Here, Espf = o p f /e p f ; a p f and s pf are respective values of fasteners at failure; F pf is the 
bolt or rivet cross-sectional area; and Cf and nf are the fastener spacing and the number of 
fastener rows, respectively. 

However, fracture may be because of bolt head breakage and/or skin tear at the bolt head; 
in these situations the c p f coefficient can be computed on the basis of data from finite-element 

analyses, similar to those in Figs 2.8 and 2.1 1 . 

With (1.29) derived, we have 

y f =|aA u (N5') 2 X 

Upon substitution of (2.30), (2.32), and (2.34), the relation (2.36) between the failure load q 


and the link fracture energy X appears: 

q p =Jx/c, 


(2.36) 


in which 


- - - - 2c pf 

C = C p f + Cbf , C p f = ~ , 


c bf = 


2c bf 
aAn ’ 


(2.37) 


If one knows parameters of the rib and fasteners (that is, 
can be utilized to calculate the relative per-unit-length load q p 


c is known) the formula (2.36) 
applied by the panel to the rib 


through the fasteners. 

In accordance with (1.28) and (1.44), the panel coming from state I into state S (with S 
being either II or HI) is loaded with the compressive stress resultant 

N = Ncr+Vf^> 

from here and (2.36), the expression is, 
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Figure 2.12 represents «-N-q p » diagrams for three versions of the stiffened panel in which the 
skin has the following characteristics: 


D? 


Dfl A 22 u\ 1 

= 2J8 ’ ^ = 3 - 64 > “^r = 4 - 63 > 


a 3 ° 


D 22 A?, ' D 3 ° - ’ Af, 

the three panels have identical values of bending stiffness ratio D n /D 2 2 = 70 and coefficients 
a 0 = 0.3 and c = 1 50 ; whereas, the panels differ in the ratio A 22 /A j j : 

A 22 


L 11 


■ = 3.44, 4.30, 5.72. 


Diagrams in Fig. 2.12 suggest that if the longitudinal stiffness of the panel stringers is increased, 
then a particular N value (that is a minimum relative load which separates the panel from ribs) 
takes place at a lower load q applied to the panel/rib fasteners. 

Upon computing the available q p values the plots in Fig. 2.12 can be employed to 
determine whether the panel separates from ribs under the compressive load. If q <(q ) 

then the joint will become broken at a load N < N™ . However, if q > (q ) then the panel 
buckles between ribs at N = N ^ . 

If one would specify rib and fastener parameters which ensure the panel to buckle 
between ribs, then the parameter c should be determined from the condition q >(q ) .Let 

us consider now the dependence of (q J on some of the panel parameters. 

Compute X = X CT by using (1.47) and (1.48) and substitute the value into (2.36) to 
derive the dependence of critical (q p )^ values on a 0 for the panel concept versions involved; 

see Fig. 2.13. For each value of A 2 2 /A n within the interval (0.3<a o <0.45} the (q ) 
variable does only slightly vary with a 0 . 

A similar kind of picture for (q p j^ as a function of D n /D 22 for the above values of 

A 22 /All > a o =0.3, and c=150 may be seen in Fig. 2.14. It is clear that (q ) does not 
almost depend on Dj ] /D 22 when this ratio is above 50. This fact allows us to derive a unified 
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relation between I q_ ) and relative compliance c of the panel-to-rib joint for a 0 2:0.3 and 
V P/cr 

Di i /D 22 ^50; the relation is demonstrated in Fig. 2.15. 

There exists a really essential interval of variation of the relative panel-to-rib joint 

compliance: 1 00 < c < 400 ; here, the critical value of the relative failure load ( q ) varies from 

V P/cr 

0.04 to 0.08 - this means that the panel-to-rib joint will fail under a load exceeding 4% - 8% of the 

— IV 

critical compressive load N CT ; note also that the stiffer the rib, the higher the failure load 

I q ) and the stronger the fasteners must be. 

The function depicted in Fig. 2.15 makes it possible to specify parameters of usual ribs 
and fasteners so as to ensure the panel to buckle between ribs without breaking the joint. 
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CONCLUSION 

The NCC-1-233 programme (in Appendix 5) dated September 15, 1998, has been the 
basis for the investigation of how compression-loaded multibay panels manufactured of 
composite materials buckle upon failure of joints with supporting elements. 

The following principal results should be noticed here: 

- a new method for studying the buckling of a multibay stiffened panel made of 
composites is proposed; it allows for destruction of panel joints with transverse and 
longitudinal supports in the course of buckling; 

- the solution to the problem is derived; parametric research on failure load dependence 
on available link fracture energy and other panel parameters has been conducted; 

- notions of critical link fracture energy and link strength for a multibay panel are 
introduced; the critical values outline the domain within which the panel buckles with 
nodal lines being on supports; 

- techniques to compute the link fracture energy for some types of joint are proposed; 

- a relation of link fracture energy to the load diagram and failure load of the joint is 
established; and 

- a methodology for calculating the support/fastener parameters necessary to ensure a 
strong joint between the panel and supports is presented; usual web-like ribs and 
fasteners are shown to be necessarily designed so as to carry a load which is 4 % to 8% 
of the longitudinal compressive critical force for the panel. 

The present study method makes it also possible ’to solve the following problems of 
practical interest: how the compression-loaded composite panel buckles upon local separation of 
stiffeners from the panel and how strong the joint between the skin and stiffeners should be for 
the separation not to occur before the general buckling takes place. 

By proceeding in a similar way, the problem of compression-loaded sandwich panel 
buckling with separation of the composite skin from the core over a certain (local) portion can 
be solved. . 

Efforts including solution of these and some related problems may be a continuation of 
the composites study programme within cooperation between NASA and TsAGI. 
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Fig. 1.4 



Fig. 1.5 
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Fig. 1.8 
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Fig. 1.11 
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Fig. 1.15 



















Fig. 1.19 
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